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Abstract 



There is a widespread belief that the general two-Higgs-doublet model 
(G2HDM) behaves unnaturally with respect to evolution of the flavor- 
changing neutral Yukawa coupling parameters (FCNYCP's) - i.e., that the 
latter, although being suppressed at low energies of probes, in general in- 
crease by a large factor as the energy of probes increases. We investigate this, 
by evolving Yukawa parameters by one-loop renormalization group equations 
and neglecting contributions of the first quark generation. For patterns of 
FCNYCP suppression at low energies suggested by existing quark mass hi- 
erarchies, FCNYCP's remain remarkably stable (suppressed) up to energies 
very close to the Landau pole. This indicates that G2HDM preserves FC- 
NYCP suppression, for reasonably chosen patterns of that suppression at low 
energies. 
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T.D. Lee proposed already in 1973 the general model with two Higgs doublets 
[G2HDM, referred often also as 2HDM(III)]. This is the model with the most general Yukawa 
couplings of the two Higgs doublets (to quarks). He investigated the model with emphasis 
on CP-violating phenomena, and was followed by Sikivie et al. [|J. 

Glashow and Weinberg |J pointed out in 1977 that only those models with two Higgs 
doublets which possess specific discrete [or U(l)-type] family symmetries in the Yukawa 
sector [2HDM(II) and 2HDM(I)] have the property of complete suppression of all flavor- 
changing neutral Yukawa coupling parameters (FCNYCP's), and this is then true at any 
evolution energy^ Since then, a large part of physics community has apparently regarded 
such models as being the only ones with a natural suppression of FCNYCP's. It appears 
that G2HDM's were then not investigated until late eighties. 

During the last twelve years, there has been a moderate resurgence of works on G2HDM's 
B~ 0-0 These works investigate low energy phenomena of G2HDM's, allowing nonzero but 
reasonably suppressed FCNYCP's at low (electroweak) energies. Cheng, Sher and Yuan 
(CSY) H and Antaramian, Hall and Rasin (AHR) |J proposed specific ansatze for these 
parameters at low (electroweak) energies, motivated largely by the existing mass hierarchies 
of quarks. They pointed out that neutral scalars can then still be reasonably light (M ~ 10 2 
GeV) without violating available data on flavor-changing neutral processes. Thus, these 
two groups showed that FCNC suppression in G2HDM's is not intrinsically unnatural (i.e., 
ad hoc), since it can be motivated by existing quark mass hierarchies. Certainly, these 
arguments regard only the structure of the theory at low (electroweak) energies of probes. 
We note that most of the phenomenological investigations in G2HDM's by other authors 
used essentially ansatze of these two groups. 

We are not aware of any work on the second aspect of G2HDM's, i.e., the question 
whether these models really behave unnaturally with respect to evolution of FCNYCP's. 
The present work attempts to address this question, by investigating one-loop evolution of 
these parameters, for the cases when the patterns of their values at low energies are described 
by CSY-type of ansatze and modifications thereof. 

Yukawa interactions of quarks in G2HDM's in any SU(2)L-basis are 

+U$ ) C$* {1) )*$ + ugHW*®)*® + h-c.} • (1) 

Parameters and quark fields are in an arbitrary SU(2)L-basis, as indicated by tildes above 
them. Superscript (E) at the Lagrangian density means that the theory has a finite effective 



1 Low energy experiments show that those flavor-changing neutral coupling parameters (FCNCP's) 
which don't involve t quark are suppressed in nature at low energies E~m q . For FCNCP's involving 
t, experimental evidence is not yet available. 

2 Refs. and Refs. Jjfl concentrate largely on CP-violating and FCNC-violating phenomena, 
respectively. 
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energy cutoff E (energy of probes). This superscript is omitted at the fields and at the 
parameters for simpler notation.f] The following notations are used: 



V2 



d k ) i -A k ) 



/0\ ^ e^/0 



Here, t> [= v(E)} is the usual vacuum expectation value (VEV) needed for the electroweak 
symmetry breaking: t>(.E ew )^246 GeV. Phase difference r] = r]2 — i]i between the two VEV's 
in (H) may be nonzero and represents CP violation originating from the purely scalar sector. 

Popular 2HDM(I) and 2HDM(II) models proposed by Glashow and Weinberg || are 
special cases (subsets) of this framework: U {1) = = in 2HDM(I), and U {1) = L> (2) = 
in 2HDM(II). In these models, mass matrices M<& an d in the mass basis are propor- 

tional to the relevant nonzero Yukawa matrices U, D, respectively. Therefore, there are no 
FCNYCP's in the physical (quark mass) basis, since off-diagonal elements are zero there by 
definition. Moreover, this remains true even when radiative corrections are included, i.e., 
when the model is evolved from a high "bare" energy E = A to low energies E~m q . This 
is so because the mentioned structure of the Yukawa sector is ensured by specific discrete 
symmetries. 

There are no such symmetries in the Yukawa sector of G2HDM , so a priori it is unclear 
whether suppression of FCNYCP's would persist when we increase energy of probes beyond 
electroweak scales. We carried out such an analysis, by using one-loop renormalization 
group equations (RGE's) for Yukawa coupling matrices appearing in flU). Those RGE's were 
derived by us earlier ||. 

In order to interpret more easily the RGE results at high energies of probes, as well as 
the initial conditions at low (electroweak) scales, it is convenient to introduce the following 
derived quantities: 

$'« = cos +e~ iri sin /5$ (2) , $ /(2) = - sin/3$ (1) + e~ iri cos/3$ {2) , (5) 

where: tan (3 = — =>■ cos /3 — — , sin (3 = — ; i] = m — f?i • (6) 

Vi v v 

VEV's of the redefined scalar isodoublets are: e^ 1 = (Q,v/V2), ($ /(2) )o = (0,0). 

Hence, isodoublet is responsible for the quark masses. Below we will see that 
leads to FCNYCP's. We now rewrite original Yukawa Lagrangian density (0) of G2HDM in 
terms of these fields 



3 For renormalized quantities, I?~10 2 GeV. 
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- E {A,(gf $ /(2) )^ } + ^(gf $' (2) )^' } + h.c.} . (7) 

G (c/) and G*^) are rescaled mass matrices, and U and Z) the corresponding Yukawa matrices 
"orthonormal" to and G( D ) 

= cos /?X (1) + e Ti '" sin /3X (2) = v^M w /« (X = U,D) , (8) 
X = - sin /3X (1) + e Tir? cos /3X (2) (X = U, D) . (9) 

Minus sign in exponents is for X = U, and plus for X = D. Transition to the quark mass 
basis (at a given energy E) is implemented by biunitary transformations involving unitary 
matrices V^ R and V^ R 



U^yVfjym. G V) = v u & u) v u^ G V) = 5 .. m M^2 /v , (10) 

D = V°DV^; GW = V?GWvg*, 4 D) = 6 ijm ^V2/v ; (11) 
u L = V?u L , u R = V R J u R , d L = V°d L , d R = V°d R . (12) 



Parameters and fields without tildes are in the mass basis. Superscripts (E) are omitted for 
simpler notation. CKM mixing matrix is V = V^V^ . Neutral part of Lagrangian density 
(0) in the quark mass basis is then 

+ D l3 dfd^f ) +^f ) ] + U^uWP-^'P] + h-c} , (13) 

where summation is performed over repeated flavor indices. We see from (|13|) that the U 
and D matrices, defined by ©-(11) through the original matrices U^> and D^> of G2HDM 
([I]), allow the model to possess in general FCNYCP's. This is so because in the quark mass 
basis only the (rescaled) quark mass matrices G^ and G^ of ( [T0|) -(|IlD [cf. also @] are 
diagonal, while additional matrices U and D are in general not. 

As mentioned earlier, CSY ]5] argued that U and D matrices in the quark mass basis 
and at low energies E have the form 

V2 r M M n I m\ /-(d) 



with: CS^QT ~ 1 for E ~ £ cw (~ M z ) . (15) 

This form, at least for the diagonal elements, is suggested by the existing quark mass hier- 
archies and the requirement that there be no fine-tuning in which large Yukawa terms 

(and: D^l) would add together via © to result in much smaller terms Ujk {Djk) - for this, 
inspect Eqs. ©— ©, but this time in the quark mass basis (no tildes). Moreover, this ansatz 
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turns out to be phenomenologically viable, i.e., the known suppression of flavor-changing 
neutral processes at low energies (not involving on-shell top quarks) can be reproduced. 
Similar (but not identical) ansatze have been proposed by the authors of pi (AHR), mo- 
tivated by their requirement that the Yukawa interactions have approximate U(l) flavor 
symmetries. 

Now we are prepared to present some typical examples of the RGE evolution of Yukawa 
parameters in G2HDM's. For simplicity of analysis, we neglect contributions of the first 
quark generation. Further, we assume that all original four Yukawa matrices D^> are 
real and the VEV phase difference r\ is zero (no CP violation). 

For the boundary conditions to the RGE's, at the evolution energy E = M z , we first took 



the CSY ansatz ©-(PD , with ^ = 1 = ^ or £g° =2 = £f , for all i,j = 1,2. We emphasize 
that 2 = 1 refers now to the second quark family (c,s), and % — 2 to the third family (t,b). 
For the (2x2) orthonormal CKM mixing matrix V we take V 12 (M Z ) = 0.045 = -V 2 i{M z ). 
Values of other parameters at E = M z were chosen to be: 

tan p = 1.0; v = <Jv\ + v\ = 246.22 GeV; a 3 = 0.118, a 2 = 0.332, a x = 0.101; m c = 0.77 GeV, 
m s = 0.11 GeV, m& = 3.2 GeV, and m t = 171.5 GeV. These quark mass values correspond 
to: m c (m c ) w 1.3 GeV, m,(lGeV) w 0.2 GeV, m b (m b ) ^4.3 GeV, and mf hys ' w 174 GeV 
[rat(rat) ~ 166 GeV]. For a^(E) we used two-loop evolution formulas, with threshold effect 
at E^m^ hys ' taken into account; for otj{E) ( j = 1,2) we used one-loop evolution formulas. 

This simplified framework resulted in 18 coupled RGE's [for 18 real parameters: v 2 , 
tan/3, Uij, £>ij, G\J\ G+f ], with the mentioned boundary conditions at E = M z . The 
RGE system was solved numerically, using Runge-Kutta subroutines with adaptive stepsize 
control (given in 0). 

The results for the ratios of FCNYCP's X^E)/ X^Mz) {X = U,D; i ^ j) are given 



for the case of ^ = 1 



Ad) 



in Fig. 1. The Figure shows that all the FCNYCP's are re- 



5.5 - 

5 
4.5 

4 
.1.5 

o 

'$ 3 
2.5 

2 - 

1.5 

1 

0.5 



(o L t R ) U 12 (E)/U 12 (M Z ) 

<s L b R ) D 12 (E)/D I2 (M Z ) 

(t L c R ) U 21 (E)/U 21 (M Z ) 

<b L s R ) D 21 (E)/D,!(M Z ) 



log 10 (Energy/GeV) 



FIG. 1. Evolution of FCNYCP ratios U ij {E)/U ij (M z ), D ij {E)/D ij (M z ) {i^j), when Qj>=C 
(for all i,j = 1,2). 



markably stable as the evolution energy increases. Even the up-type FCNYCP's, appearing 
in couplings involving the heavy t quark, remain rather stable. Only very close to the t- 
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quark-dominated Landau pole (~10 13 GeV) these parameters start increasing substantially. 
For example, in the down- type FCN sector (b-c) the corresponding ratio D 2 x{E) / D 2 \{Mz) 
acquires its double initial value (i.e., value 2) at E ~ 0.7E po \ e , which is very near the pole. 
About the same is true also for U 2 i(E) /U 2 i(M z ) ■ For the ratio Di 2 (E) / D 12 (M Z ) the corre- 
sponding energy is even closer to £^ po i e , while for for the t-quark-dominated U\ 2 {E) J U\ 2 {M Z ) 
it is somewhat lower. We draw the conclusion that the case ^ = 1 = ^ shows no "un- 
naturality" concerning the behavior of FCNYCP's at high energies. We may also compare 
this stability with that of other ratios. For example, in Ref. || we showed also behavior of 
those ratios of neutral Yukawa parameters which don't involve flavor-changing couplings: 



Ujj(E)/Ujj(M z ) and DJ E) : I) u iM x ): G%\e)/G%\M z ) and G%\E)/G%\M Z ). We 
saw that several (but not all) of these coupling parameters are also reasonably stable. The 
results for FCNYCP's of Fig. 1 are interesting and perhaps surprising. We should bear in 
mind that the off-diagonal Yukawa parameters appearing in Fig. 1 are at low energies by 
several factors smaller than the third generation diagonal parameters, due to CSY ansatzfj 
Therefore, the fear that the latter, substantially larger, parameters would "pull up" the sup- 
pressed off-diagonal ones (by a large factor, or even by orders of magnitude) as the energy 
of probes increases, is intuitively justified. Fig. 1 says that this doesn't happen. 

One may object that behavior shown in Fig. 1 is due to a special choice C,^ = 1 = £,if- 
In fact, in Ref. 13, where the relevant RGE's were derived, we performed numerical analysis 



only for the case Qj 



(a) 



Ad) 



We have recently performed calculations for variations 



of the CSY ansatz. Results are independent of the VEV ratio tan j3 (cf. also discussion 



in P|). In Fig. 2 we show results when £, 



(«) 



Ad) 



at E = M z (for all i,j = 1,2). The 
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FIG. 2. Same as in Fig. 1, but for the choice = 2 (for all i,j = l,2). 

Landau pole is now of course lower (~10 TeV), but the behavior with respect to the stability 
of FCNYCP's remains qualitatively the same. When we vary some of the £y parameters, 



Note that D\ 2 = D 2 \ at E = M Z is two orders of magnitude smaller than U 22 . 
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the results remain qualitatively the same, while the Landau pole is influenced largely by 
the (t-quark-dominated) up-type CSY parameter £22 [Z7 2 2 {Mz), cf. Eqs. ([14]) -(115)1. We 
investigated also cases which go beyond CSY ansatz. For example, we chose to suppress the 
up-type off-diagonal parameter even more drastically, by taking ^ = £21^ = 0.05163 and all 
other £y parameters equal to 1. For such a choice, we have Di 2 = D 2 % — U\i — U21 at E = Mz- 
The results are given in Fig. 3, and are very close to those of Fig. 1. 
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FIG. 3. Same as in Fig. 1, but for the choice =0.05163 (other fa's are 1). 

We conclude that in G2HDM's, with effects of the first light family neglected, the flavor- 
changing neutral Yukawa coupling parameters (FCNYCP's) remain remarkably stable as the 
energy of probes increases up to the vicinity of the (top-quark-dominated) Landau pole. 
The conclusion is insensitive to any specific variation of the CSY ansatz flnp-flTSD for Yukawa 
parameters at low energies, and apparently survives even when going beyond this ansatz by 
suppressing low energy FCNYCP's even more. Thus, in the presented framework, the fear 
that the suppressed FCNYCP's at low energies are pulled up drastically by the much larger 
diagonal Yukawa coupling parameters as the energy of probes increases doesn't materialize. 
The described framework appears to behave naturally in this respect, FCNYCP's remain 
remarkably suppressed even at higher energies although there is no explicit exact symmetry 
which would ensure complete suppression of these parameters. In this sense, we have an 
indication that the G2HDM's are not unnatural, in contrast with the widely held beliefs. 

One may still argue that our conclusions might change when contributions of the first 
quark generation are included. We intend to perform this extension of numerical analysis in 
the near future. 



Abbreviations frequently used in the article: 

CSY - Cheng, Sher and Yuan; FCNC - flavor- changing neutral currents; FCNYCP - flavor- 
changing neutral Yukawa coupling parameter; G2HDM - general two-Higgs-doublet (Stan- 
dard) Model; RGE - renormalization group equation; VEV - vacuum expectation value. 
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